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ABSTRACT 


This thesis examines the representation of dimensional units as prime numbers 
to perform dimensional analysis within a computer-based model management system. 
A computer program applies this concept to simple span structural steel beam design, 
an engineering stress and strain problem. Most common applications of computers 
manipulate only the numeric value of the measure of physical objects. The user 
manually ensures that data is processed according to the meaning of its units. Prime- 
encoding of dimensional units in this application provides a numeric method of 
validating dimensional consistency in mathematical expressions for use on a 
computer. This study is implemented in TEFA, a computer-based modeling system with 
an embedded Prolog programming language. The beam design application 
demonstrates that model representation using prime-encoding of dimensional units 
simplifies the overhead required in data manipulation, and helps maintain 


meaningful results in the numerical processing of data. 
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I. INTRODUCTION 


A. OBJECTIVES 

This thesis examines the representation of dimensional units as prime numbers 
to perform dimensional analysis within a computer-based model management system. 
A computer program applies this concept to simple span structural steel beam design, 
an engineering stress and strain problem. Most common applications of computers 
manipulate only the numeric value of the measure of physical objects [1: p. 478]. The 
user manually ensures that data is processed according to the meaning of its units. 
Prime-encoding of dimensional units in this application provides a numeric method of 
validating dimensional consistency in mathematical] expressions for use on a 
computer [2: pp. 2-3]. This approach has been implemented in TEFA, a computer -based 
modeling system with an embedded Prolog programming language [3]. The beam design 


application demonstrates that: 


Model representation using prime-encoding of dimensional units sümplifies the 
overhead required in data manipulation, and helps maintain meaningful results 
in the numerical processing of data. 


Bradley [4: pp. 403-404] suggests that little attention has been given to 
representing and verifying data types used by computer programs. In most modeling 
applications, a computer-based mathematical model of a physical process typically 
manipulates the numerical portion of the model. However, the program may not be 
capable of distinguishing qualitative differences between the data it manipulates, 
making the program prone to errors. The goal of the modeler is to process numerical 
data according to its meaning, including constraints placed upon the data by the model. 
Bradley proposes that a system consisting of descriptions of model variables and rules 


for manipulating them can address this problem. 


A simplified problem of structural steel beam design is used to demonstrate the 
representation of dimensional units as prune numbers in dimensional analysis. Beam 
design. in this thesis, is represented as a series of mathematical models describing the 
physical phenomena of a load upon a beam. In this context. a mathematical model is 


specified by inputs. objective functions, and constraints [5: pp. 2-7]. 


B. BENDING OF DUCTILE STEEL BEAMS 

Parker [6: p. v] proposed that there is a continuing need to simplify structural 
steel design procedures for some types of butldings. Architects and engineers, 
henceforth referred to as users. require knowledge of engineering equations used for 
structural steel beam design to effectively oordinate design and constm ‘ion. This 
includes validation of data available to the user for use in the appropriati *quation. 
The beam design example is used to examine dimensional validation typically 
performed in an engineering problem. 

Beam design begins with an examination of spanning and loading conditions of a 
building structure. Parker [6: p. 45-48] states that loads supported by beams are 
classified as either concentrated or distributedloads. A concentrated load is assumed to 
act at a point since it extends over a relatively small a portion of the beam length as 


shown in Figure 1. 


Figure 1 Simple Beam with Concentrated Load 


A distributed load on a beam extends over a substantial portion of the entire 
length of a beam. Distributed loads are commonly referred to as uniformly distributed. 
The load has a uniform magnitude for each unit of length, such as pounds per linear 
foot or kips (1000 pounds) per linear foot, as shown in Figure 2. 


Figure 2 Simple Beam with Uniformly Distributed Load 


The focus of this application is the design of a simple beam with a uniformly 
distributed load. A general procedure for this case of beam design is briefly outlined as 
follows [6: pp. 130-132]: 


Step 1: Compute the load, W , in kips, that the beam will be required to support. If w 
is the uniformly distributed load in kip/feet, and L is the length in feet, then 


W=wL. 
Step 2: Compute the reaction, R, in kips: 
R=W / 2. 
Step 3: Compute the maximum bending moment, M, in kip-feet: 
M=WL/8. 


Step 4: Determine section modulus in inches? by use of the beam formula where 
F, = 24 ksi, the allowable bending stress of A36 steel: 


Step 5: Refer to the table of structural shapes, Properties for Designing, and select a 
beam with a section modulus S equal to or greater than that required by Step 
4. See Appendix A. The compact section criteria in the table should confirm 
that b, / 2t,, is equal to or less than 8.70, the limiting value for A36 steel 
where F,=24 ksi. The beam selected will be adequate for bending stresses 
subject to checking of section classification (compact or noncompact) and L, 
or L, controls. 


Step 6: Check the beam for shear by confirming that the computed shearing stress 
f,= V/A, in ksi is less than the allowable shear stress F,,= 0.40 x 36 = 14.5 


ksi for A36 steel. The static vertical shear on the beam, V in kips, is equal to 
R. The area of the girder web is determined by A, =dt,,. in square inches, 
where d is the depth of the beam in inches and t, is the web thickness of the 
beam in inches. 


Step 7: Verify that the maximum actual deflection caused by the loading is less than 
the allowable deflection. Actual deflection is defined as follows: 


A=(5xfx12)/(48xExc). 


This simp.ties to: 
A 902353 L2)/ a: 


where A is in inches, f is in ksi, 1 = L x (12 in/ft) is in inches, E is the modulus 
of elasticity of steel equal to 29,000 ksi, and c is the distance from the 
neutral axis to the extreme fiber, where c = d/2 in inches. Allowable 
deflection A. in inches. is equal to 1/360 of the span or (L x 12)/360. If the 
actual A exceeds the allowable A, deflection may be solved for I, moment of 
inertia where I is in inches, required to limit A to the allowable value. This 
part of the procedure will not be attempted in this study. 


Appendix B provides an example of a beam design problem using the above 
procedure. Appendix C provides a reference of general nomenclature used by the 
problem, and Appendix D provides a list of the abbreviations referenced in the design 
procedure. The purpose of the beam design application is to verify the dimensional 
consistency of the required mathematical express ns. The application also solves 
the expressions, recommending a wide flanged steel beam with a physical 


configuration fulfilling the design criteria of the problem. 


C. DIMENSIONAL ANALYSIS OF BEAM DESIGN 

Beam design illustrates typical difficulties of analyzing engineering expressions. 
Unit conversions are often required, as in the case of determining section modulus. The 
computed value for S, using the beam formula Ss M/F,. must be in inches? for 
comparison with section modulus values of the beam table. If the value of M is entered 
with dimensional units of kip-feet then M must be converted from kip-feet to kip- 
inches to usefully represent S in inches?. The dimensional units of F, can be 
represented as either ksi, kip/inches?, or 1,000 pounds/inches?, depending upon how 
data is presented to the user, or upon the user's preference. With these problems in mind, 
a fast, accurate, and automated system which performs checks of dimensional 


consistency is desirable. 


D. OUTLINE OF THIS PAPER 

In Chapter II, the concepts of measurement, physical algebra, dimensional 
analysis and application to modeling languages are briefly be discussed. Chapter III 
examines laws of dimensional consistency necessary to perform dimensional analysis 
of mathematical expressions. Prime-encoding of dimensional units is introduced for 
the purpose of performing dimensional analysis. Chapter IV presents the application of 


dimensional analysis and prime-encoding of dimensional units to the steel beam 
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design problem within a modeling language. Chapter V summarizes the results of the 
study and presents recommendations for further study. Appendices A through D 
provide a sample beam design problem with associated beam tables, nomenclature and 
abbreviations. Appendices E and F contain prime-encoded dimensional units and rules 
for determining dimensional consistency. Appendices G through J contain the beam 
design application, support processes and sample output. Appendix K describes the 


computer platform used to develop the application. 


II. BACKGROUND 


A. MEASUREMENT AND PHYSICAL ALGEBRA 

Massey [7: p. 11-12] states that measurement is basically a comparison of things 
of the same kind. Measurement uses a unit defined as a standard amount of a quantity 
to compare with another of the same kind of quantity. The comparison is a magnitude 
with a numeric and a unit component. For example, to express the length of an object, a 
number quantifies the number: tandard units equal to the length of the object, and a 
unit identifies the standard amc it of that quantity. 

The algebra used to describe describe physical objects, as suggested by Massey [7: 
pp. 12-13], is different from the algebra of pure mathematics. The algebra of pure 
mathematics, or ordi: ary algebra, is the expression of relations among numbers, 
though symbols are often used to represent numbers. Physical algebra describes 
relations among the magnitudes of physical quantities such as force, velocity, mass, 
energy, and so on. Physical algebra describes how the magnitude of one quantity 
depends on the magnitude of the others. Mathematical operations of addition, 
subtraction and comparison are restricted to quantities of the same kind. For example, 
a mass and an interval of time cannot be compared, but one mass and another mass 
can. 

Masscy [7: p. 13] states that physical algebra expresses relations between 
magnitudes of the same kinds of physical quantities. Magnitudes of the same kind of 
physical quantities are also expressed in terms of the same dimensional unit. When 
both sides of an equation have the same units, the equation is dimensionally 
homogeneous. 

Relations in physical algebra, Massey [7: p. 13] asserts, must have two kinds of 
consistency. First, there is the quantitative relation of ordinary algebra in which both 
sides of an equation are compared in numerical magnitude. Second, there is a 
qualitative relation such that terms that are added, subtracted, or compared represent 
th: ime kind of quantity. Ordinary algebra is basically a means of comparing 
nu. rs. Aerea. ^hysical algebra is a means of comparing the magnitu es of similar 


qu des. 


B. DIMENSIONAL ANALYSIS 

Langhaar [8: p. 1-4] describes dimensional analysis as a method of deducing 
information about physical phenomena if it can be described by a dimensionally 
correct equation. Scientific reasoning of the physical world is based on concepts of 
various abstract entities, such as force, mass, length, time, and so on. Each of these 
entities, or dimensions, may be assigned a unit measurement and is considered to be 
independent of the others. 

Bhaskar [9: p. 73-74] states that dimensional analysis has been used in 
engineering for purposes of modeling and similitude. With dimensional analysis, 
reasoning about a system is possible without explicit knowledge of the physical laws 
that govern it. The model of the system requires knowledge of only relevant physical 
variables and their dimensional representation. Dimensional representation of 
physical variables contains a significant amount of knowledge. 

Physical representations of variables, as stated by Bhaskar [9: p. 73-74] have both 
numerical and symbolic components. The numerical component is the value of a 
variable measured in a system of units. Reasoning about the numerical component 
represented by a physical variable is constrained by the physical context the value may 
take. The symbolic component in qualitative physics is the dimensional 
representation of the physical variables. For example, in the dimensional notation of 
physics, force is usually represented as pounds x feet x seconds ^. Dimensional 
representation of a variable is also constrained by a set of laws. Dimensional 
homogeneity is the most familiar of these. One of the most widely used results of this 
concept in dimensional analysis is Buckingham's []-Theorem, proved by Buckingham 
[10: pp. 345-376] in 1914. This theorem is used to establish the number of independent 
dimensionless numbers required to describe a given physical context 
[11: pp. 3918-3919]. Langhaar [8: p. 18] explains that according to Buckingham's [I- 
Theorem, an equation that relates dimensionless products is dimensionally 
homogeneous. Stated another way, the form of the equation is independent of the 


fundamental units of measurement. 


E DIMENSIONAL ANALYSIS USING A COMPUTER 
Various computer -algebra systems and executable modeling languages [1: p. 478]. 
[3]. [4]. [11]. [12]. [13] have been proposed to include symbolic physical units in computer 


calculations. These techniqu offer automatic <“ tion of dimensionally 
inhomogeneous formulas and .tomatic conversion o! .nconsistent units in a 
dimensionally homogeneous formula. 

Stoutemyer [1: p. 478] states that attempts to include physical units in 
mathematical expressions are sometimes abandoned since conventional programming 
languages generally deal with pure numbers rather than physical quantities. The user 
of the program must manually check data for its dimensional type to ensure that the 
results of the computation are valid. Errors of dimensional consistency are difficult to 
prevent and detect in traditional programming languages. 

Stoutemyer [1: p. 478] suggested that since many programming languages have the 
ability to declare the precision of numerical variables and whether it is fixed or 
floating point, it should be possible to extend this technique to the declaration of 
variables such as units; i.e., meter, mile/hour, dollars, etc. A modeling language 
translator then checks expressions and assignments to variables for dimensional 


consistency relieving the user of the tedious, error-prone unit conversion process. 


D.  ALTERNATE METHODS 

There are several computer programs incorporating dimensional analysis. 
Though they share the similar goal of validating dimensional consistency of a 
mathematical expression, they are quite different in their approach. The survey of 
alternative methods presented below is not exhaustive, but illustrates three different 
techniques of representing ar extended data type  »ntaining numerical and 
dimensional information. 

The first method, as described by Barnes [14: p. 3-14], 1s a database representation 
using a semantic type checking system for use with relational databases. The unit 
component of a variable, or semantic information, is associated with its numerical 
component in a data dictionary. This extended numeric data type, defined by concept, 
quantity and dimension, consists of a value description and a semantic description. 
The value description of the data type consists of a number and its unit of measure. The 
semantic description consists of a quantity and a concept, which is an object attribute 
possessing that quantity. Database queries are verified for dimensional consistency 
and dimensional units are converted across systems of measurement as necessary. 


Concepts are used to construct a concept hierarchy ensuring that queries are consistent 


with the semantics of the database. This abstract numeric data type is used with the 
relational data model for dimensional validation and unit conversion. 

Another variant of a computer-based dimensional analysis system proposed by 
Hirschberg [15: p. 2-9] uses the Buckingham Pi-Theorem [10] to order sets of equations 
by their importance. This system is implemented in SYMBOLANG, a list structured 
symbol manipulator based in the FORTRAN programming language. The symbol 
manipulator operates on a set of symbols rather than numbers. For example. a symbol 
manipulator multiplies N+1 by N-1 to obtain N?-1. A list structure stores and 
manipulates data by defined relationships. Solution of the Pi-Theorem involves the 
formation of Pi-Terms ordering the equations by their importance. Once a Pi-Term is 
formed, the dimensionless set of terms is solved as a linear set of equations. The 
numeric solutions are then paired with their associated parameters. Hirschberg 
proposed that an extension of this work would include checking Pi-Term solutions 
against a set of well known dimensionless numbers. Data for the parameters would 
serve as input producing numeric answers. Hirschberg considers the speed of the 
computer an important advantage because it generates many permutations of 
equations yielding a large number of solutions, which can be solved by regression 
analysis for a best fit solution. 

Stoutemyer [1: pp. 479-480] proposed yet another method using a computer -based 
symbolic algebra system to process dimensional units. The computer algebra technique 
is implemented in MACSYMA, developed by the MATHLAB Group. The method extends 
a programming language's ability to declare not only numerical variables, but also 
variables with a dimensional unit component such as meters, miles/hour, dollars, etc. 
This technique detects dimensionally inconsistent formulas and converts inconsistent 
units within a dimensionally homogeneous formula. A translator checks expressions 


and variables assignments for dimensional consistency. 


E. DIMENSIONAL INFORMATION IN MODELING LANGUAGES 

Having identified the need for verifying dimensional consistency of an equation 
describing physical relationships, it would be useful to have a flexible, easy to develop 
method of representing mathematical expressions describing the physical situations. 
Fourer suggests (13: pp. 143-155] that algorithmic solutions to solve problems within 


computer programs are explicit rather than symbolic. They are designed for 


convenience and efficiency of handling data in a computer program rather than for 
clarity to the modeler. As a result of this approach, verification and modification of the 
modeler’s intentions become a problem of debugging computer programs. In addition, 
the description of the problem is highly dependent upon the form of the algorithm 
chosen to represent it. 

Fourer [13: pp. 155-163] has argued modeling languages offer several advantages 
over algorithmic descriptions of a problem. The modeling language is not a 
programming language, but a declarative language that describes the modeler's 
intentions in a form that can be interpreted by a computer. The ability to describe the 
problem in the form of its mathematical model provides inc -ndence of any 
particular algorithmic form. Verific: a is reduced to the tas debugging the 
modeling language's representation o. he problem. Bradley [12: p. 27] states that 
present research focuses on the use of modeling languages as executable © npute” 
programs. In general, modeling languages can be used to implement type caiculus ius 
dimensional systems. A variable in the model is assigned a type that consists of its 
concepts, quantities, and units of measurement. This representation permits checking 
the composition of expressions, and defining a hierarchy of concepts with inheritance 
of properties. 
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III. DIMENSIONAL ANALYSIS 


A. INTRODUCTION 

Dimensional analysis, as asserted by Massey [7: p. 108] and Langhaar [8: p. 1], 
cannot by itself provide a complete solution to a problem, but can provide a means for 
simplifying complex problems. Laws of dimensional consistency provide a means to 
check the validity of mathematical expression. Dimensional arithmetic provides a way 
to manipulate dimensional units. Together with prime encoding of dimensional units 
terms, as described by Bhargava [2] it is possible to develop a numeric process for 
executing these techniques in a computer program. 

It is shown by Bhargava [2] that dimensional manipulation can be viewed as 
numerical arithmetic by recognizing the nature of dimensional arithmetic. Prime 
encoding supports this by representing each fundamental dimensional unit as a prime 
number. By the unique factorization theorem from number theory, numeric arithmetic 


applied to the prime encoding system follows the laws of dimensional arithmetic. 


B. DIMENSIONAL ARITHMETIC OF UNIT MEASURES 

Dimensional arithmetic provides a technique for manipulating the dimensional 
component of information [2: p. 3]. Dimensional arithmetic operates on dimensions in 
a similar manner to the way arithmetic operates on numbers. 

A variable in a mathematical model may be defined by the quantity it measures 
and its dimensional unit. For example, in the beam problem, the section modulus of a 
beam, S in inches?, is equal to moment, M in pound-inches, divided by bending stress of 
steel, F, in pound/inches?: 

S=M/F,. 


The quantity of the variable generally has a base unit of measurement and may 
sometimes contain other dimensional units [4: p. 4]. These units are related to each 
other by laws of conversion within a system of measurement, such as one foot is equal 
to 12 inches. Dimensionless quantities without units, such as constants or ratios, may 


be represented simply by the number 1. 
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The International Metric System, SI, illustrates a standar1 system of seven 


fundamental quantities and base dimensional units [16]. 


UNITS FOR THE SI SYSTEM OF MEASURES 


Quantity Measured Base Unit Other Units 

Length Meter Kilometer, Centimeter 
Mass Kilogram Gram, Milligram 
Time Second Hour, Minute 
Ampere Electric Current 

Temperature Kelvin Celsius, Fahrenheit 
Luminous Intensity Candela 

Amount of Substance Mole 


Dimensional units [2: p. 5-6] may be classified as fundamental or derived units. A 
unit is considered to be a fundamental unit if it is not a ; ict of oti: . units. A derived 
unit is product of other units. 

Validation of dimensional information, as stated by Bhargava [2: pp. 3-6]. requires 
dimensional manipulation which supports dimensional simplification and 
verification of the dimensional equivalency of expressions. Dimensional 
simplification reduces a mathematical expression to its simplest form. For example, 
the expression for section modulus, (pound x inch)/(pound x inch?), simplifies to inch?. 
Verification of dimensional equivalency requires recognition that one expression is 
dimensionally equivalent to another. For example, (pound x inch)/(pound / inch?) is 
dimensionally equivalent to pound x inch x (1/pound) x (1/ inch) x (1/inch). 

In summary, t! following laws of dimensional consistency can be used for 


dimensional validation of expressions [2: pp. 6-7]: 


1. Two expressions may be added or subtracted only if they are dimensionally 
equivalent. 


2. Two expressions may be compared for equality or inequality only if their 
dimensions are equivalent. 


3. . ¥O expressions may be multiplied regardless of their dimensions. 

4. An expression can be inverted regardless of its dimension. 

5. The exponent of an expression must be dimensionless. It may be a fraction only 
if each dimensional unit in the expression has a power that is a multiple of that 
fraction, or if the expression is dimensionless. 


6. Transcendental functions can be applied only to dimensionless expressions. 
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Dimensional arithmetic [2: p. 7] is used to perform the Manipulations described 
above. The dimension of the sum or difference of two expressions is the same 
dimension of either only if the expressions have equivalent dimensions. If the 
expressions have different dimensions, the result is not defined. The dimension of the 
product or quotient of two expressions is the product or quotient of the dimensions of 
the two expressions respectively. An expression without a dimension has a dimension 


equal to 1. This is the dimensional multiplication identity. The dimension of an 


expression's exponent is the exponent of its dimension. 


C PRIME-ENCODING OF DIMENSIONS 

The ability to represent dimensional unit terms as prime numbers provides a 
simple way of describing mathematical expressions within a computer-based modeling 
system. Dimensional unit terms can be uniquely represented as prime numbers, and 
derived dimensional unit terms can be represented by a combination of prime numbers. 
[2: pp. 9-10] The prime encoding of dimensional units can be illustrated by arbitrarily 
assigning fundamental dimensional units in a one-to-one correspondence with prime 
numbers. Since the number 1 is used as an identity, the sequence of prime numbers for 
this example will begin with 2. The numbers 2 and 3 are assigned to dimensional units, 
pound and inch, respectively. 


Consider the expression for section modulus from the beam design problem: 
S=M/F, 
or in dimensional unit terms only, 
inches? = (pound-inches) / (pounds/ inches?) 


where S is section modulus in inches?, M is moment of the beam in pound-inches, and 
F, is bending stress in pounds/inches?. The dimensional terms of M/ F, simplify to 
inches®, which verifies dimensional consistency. By use of the unique factorization 
theorem, also known as the fundamental theorem of arithmetic, it can be shown that 
dimensional equality is equivalent to numeric equality when the unit terms are 
expressed as primes numbers. [2: p. 10-11] For example, by substituting prime numbers 
for the dimensional units of M and F, the expression for section modulus becomes: 
25—12x 97 (9/22). or B= 5. 
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D. SUMMARY 

The application of prime number encoding of dimensional units, as proposed by 
Bhargava [2: pp. 11-12], suggests a method for performing dimensional manipulation. 
Prime numbers are arbitrarily assigned to a set of fundamental dimensional unit 
terms. The prime number 1 is reserved as the identity for multiplication. Dimensional 
multiplication and division are treated as numeric integer multiplication and division. 
Dimensional addition and subtraction are simply checked for equivalency. This yields 


a relatively simple method of dimensional manipulation. In summary, dimensional 


consistency and validity of an expression may be checked by computing the 
dimensional value of an algebraic e: ression, and byc — outing the numeric value of 


the dimensional expression. 
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IV. IMPLEMENTATION 


A. MODEL DEVELOPMENT 

The beam design problem incorporates dimensional analysis through use of a 
modeling language. Modeling languages support development, documentation, and use 
of mathematical models [3: p. 1], [13: p. 144]. The modeling language chosen for this 
problem, TEFA [3], can represent mathematical models consisting of definitional 
equations. A modeling language approach permits concepts to be divided into modules 
easily shared with other applications. For example, a measurement system consisting 
of prime-encoded fundamental dimensional units, and the rules for dimensional 
consistency may be easily incorporated into other models describing other kinds of 
physical phenomena. See Appendices E and F [17]. 

Model development in a modeling language is similar to development of other 
abstract mathematical models. The purpose and the objectives of the model are clearly 
defined. Relationships between variables are stated mathematically. The 
mathematical models describing beam design can be conceptualized using a framework 


of variables and the mathematical relationships between them. 


B. BACKGROUND 

A practical outcome of the prime encoding of dimensional units, as suggested by 
Bhargava [2: pp. 2-4], [17], is a series of computer-programmable statements for testing 
dimensional consistency of mathematical expressions. Dimensional consistency is 
verified by performing the numerical equivalent of dimensional manipulation 
described by the expression. If an expression is dimensionally valid, the model 
described by the expression is dimensionally consistent. If the model is dimensionally 
consistent the test outlined in Appendix F returns a logical value of true. Otherwise, if 
the model is dimensionally inconsistent, the test returns a false value. 

Numerical operations necessary to perform a dimensional consistency test 
requires recognition of dimensional operators. It also requires association of the 
dimensional unit terms with prime-encoded equivalents. The prime-encoded terms are 
operated upon according to the rules of dimensional manipulation. To exploit the 
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prime number representation of dimensional units, a function is defined which accepts 
a dimensional expression and returns a numeric value. See Appendix G. This value is 
computed from the dimensional operation on the prime-encoded dimensional unit 


terms: 


numeric function(dimensional expression) — numeric value. 


Assuming a database of prime values for fundamental units, it is possible to 
relate the prime-encoding of dimensional units. To implement the concept of prime- 
encoding of units, a table is constructed with fundamental dimensional unit terms with 
arbitrarily assigned prime numbers. See Appendix E [17: pp. 120-126]. The prime 
number is assigned to the dimensional ur.: by defining an additional rule for the 


numeric function: 
numeric function(unit, prime number) > unit code(system, prime number, unit). 


The numeric function returns the prime number indexed by the dimensional unit term. 
For example, assume that unit code(english measurement system, 2, inch). By numeric 
function(inch, prime number) the prime number value of 2 is returned. 

This numeric function also includes rules for dimensional manipulation of the 
prime encoded dimensional units. To illustrate a general case of how this function 


performs dimensional addition: 


numeric function (expression,. ^ ), 
numeric function(expression,. :_,), 


Di ian 
numeric value is n, 
— numeric function(expression, + expression,, numeric value). 


The numeric function evaluates the expressions for the prime number equivalent of 
their dimensional units. Note that this dimensional addition rule for the numeric 
function obeys the dimensional consistency laws. In a similar fashion, the remaining 
rules for subtraction, multiplication and division of dimensional unit operations are 


developed. See Appendix G. 
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An additional tool is necessary to process dimensional expressions into a form 
that can be analyzed for consistency or validity. This tool is a list function which 
processes each term of an expression. For example, a list processor for adding 


dimensional units where “+” is the dimensional addition operator has a form: 


sum list([element, | element], answer) > 
list(element, , nj). 
list(element,, n). 
answer isn, + IM. 


The term, element,, can also represent a list. This operation continues recursively until 


all elements of the list have been processed. 


C. EXAMPLE: BENDING OF DUCTILE STEEL BEAMS 
A primary objective of beam design is to make a preliminary selection of a beam 
with the smallest section modulus larger than the section modulus computed by the 


beam formula: 
S =M/F,, such that F,= 24 kip/in? for A36 steel. 


To illustrate the implementation of this problem into a modeling language, the 


beam formula is expressed as a definitional equation: 


definitional expression(compute section modulus, S + (M * 12) / Fj). 


The variable M, computed moment of the beam, is assumed to be defined in 
dimensional units of kip-feet. The definitional expression includes the conversion of M 
into kip-inches. 

Qualitative information about the variable is described within a local variable 
declaration within the mathematical model. This declaration provides the variable 
name, the model name in which the variable is used, the local variable name within the 
model, information on whether the variable is computed internally or externally to the 
model, and also the storage type of the dimensional units of the variable. For example: 


local variable(computed section modulus, compute section modulus, S, inch’). 
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This information is mapped to a declaration to be represented by: e variable. This 
notion is referred to as quiddity, or the essential quality of the variable. In the example 
of section modulus, this information tells the modeler that moment is a derived 


dimensional unit of inch3: 


quiddity variable(computed section modulus, volume, inch?). 


Information defined in these expressions can be shared with other models in a 
knowledge base. The measurement system in Appendix E is referenced to establish 
prime encoded dimensional units. Statement of the laws for dimensional consistency 
in Appendix F is referenced to verify the dimensional consistency of the mathematical 
expression defining section modulus for a beam. This verifies that the units in the 


mathematical expression are consistent with the way they are stored in the model. See 


Appendices I and J [3]. 


D. SUMMARY 

The development of a model which incorporates dimensional validation, as 
described by Bhargava [3: pp. 7-8], consists of defining the mathematical equation, its 
variables and dimensional units, encoding the dimensional units with a set of prime 
numbers, and finally evaluating the numerical equivalent of the expression by of the 
laws definir dimensional consistency Fach equation of the beam design process is 
expressed a a mathematical m^del c« -isting of mathematically defined ^quations 
and variables. The model i- conci: :y des. ribed and documente... simplifying its use 9y 
other models. Dimensional consistency is checked by the prime-encoding of the 
definitional equation’s dimensional units. This technique verifies that variables used 
by the definitional equation are dimensionally valid. 
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V. RECOMMENDATIONS AND CONCLUSION 


A. CONSTRAINTS 

An important constraint of testing for dimensional consistency, as outlined by 
Bhargava [2: p. 15], is that the application selected for modeling must be subject to the 
standard laws of dimensional consistency. The beam design problem requires only 
simple substitution of prime-encoded dimensional units for validation of the subject 
mathematical expressions. The dimensional unit of the expressions describe physical 
measurements. A set of rules for the concept of “steel beam”, or prime-encoding of “steel 
beam” types was not considered necessary. Abstract concepts, such as those suggested by 
Bradley [12: pp. 279-280]. were not required by this problem to obtain a solution, so they 
were not explored in this study. 

In addition, Dym [18: p. 1] states that the use of a model also places another 
important constraint on the problem solving process. Models are simply a 
mathematical representation of objects and processes. A mathematical model should 
not be confused with the reality of the physical problem. Mathematics performed in the 
context of a model are an abstraction of the actual physical problem which the model 
describes. 

A constraint of the application itself is that it describes only simple span steel 
beams with a uniformly distributed load. Parker (6: pp. 45-84] proposes eight general 
cases of beam loading including concentrated loading, non-uniformly distributed 
loading, and cantilevered beams. Each case requires that different moments, shear, and 
deflections be considered for analysis. Stress parameters were considered only for A36 
steel. Though A36 steel is the most commonly used material in steel frame building 
construction, there are several other steels available, each with their own allowable 
stress values. Additionally, only a small subset of fifteen W-shapes was considered for 
testing of the model. There are over sixty W-shapes available, as well as S-shapes, M- 
shapes, channels, compound angles, and bar-joists which can serve the same general 
function of a beam. 

This system is implemented in Prolog and runs on Apple Macintosh II series of 


computers. TEFA [3]. the model management system, is implemented in Prolog. See 
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Appendix K for additional information on the hardware and software utilized in this 
study. 


B. RECOMMENDATIONS 

The steel beam design application could be extended to overcome the constraints 
of the initial model. Additional models may be incorporated into the application 
covering other beam span conditions. This system could also be applied to an 
application serving as an engineer's aid in equation checking. This technique also has 
uses in fields such as operations research and finance for analyzing and checking 
parameters of various optimization models. 

It would also be worthwhile to examine th: =asibility of using dimensional 
analysis and prime-encoding of dimensional uni: n the background of another 
application, transparent to the user. For example. 4 CAD (Computer Aided Design) 
system may contain an expert system integrated into the application, providing 
interactive equation checking and dimensional validation. 

This concept of dimensional checking may also be applied to other languages, 
such as conventional programming languages. It would be useful to examine the 
program structures necessary to parallel prime-encoding techniques of the Prolog 
language and dimensional consistency laws. This method may also be adaptable to 
spreadsheets by use of look-up tables and conditional statements describing the laws of 


dimensional consistency. 


C . CONCLUSION 

Bhargava [2: p. 14] points out that dimensional manipulation can be performed 
symbolically by existing symbolic mathematics programs. Bhargava proposes that the 
prime-encoding method of dimensional manipulation performs as well as other 
symbolic methods in retaining information about manipulated symbols. Units of 
measurement can be reconstructed from the numeric value representing the unit. This 
method is easy to implement in almost any programming language. making 
dimensional analysis possible in a wide range of applications. Since this system 1s 
designed to validate dimensional consistency of expressions, its application extends 
beyond numerical processing and database searching. The problems of dimensional 
consistency and dimensional validity are handled by the system so that the user is free 


to use a particular model across other systems. The validity of the expressions 
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contained in the models are checked, and unit systems are converted as necessary to 
permit the passing of parameters between models. This ability addresses the long range 
problems of larger databases and maintenance of consistent data. 

This thesis illustrates the feasibility of providing a computer-based numerical 
method for validating the dimensional consistency of a selected set of engineering 
equations. A fast, automatic method of validating the dimensional consistency of 
equations is demonstrated. Though one motivation of this system was to show the 
ability to convert units from one measuring system to another, it was not examined in 
this thesis. However, this system shows the usefulness of implementing a prime number 
encoding scheme for representation of units of measure in mathematical expressions 
using the inheritance properties of Prolog-type data structures. If governing laws for 
manipulation of dimensional formulas are not used, contradictions may result; for 
example, in? 2 Ib x in? [2], [19]. Such contradictions can be avoided by using standard 
conventions and techniques of analyzing dimensions recognizing the origin and 
derivation of formulas. Adherence to dimensional validation provides better 
understanding of models, and origins of mathematical expressions, and improvement 


of model representation. 
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APPENDIX A 


PROPERTIES FOR DESIGNING 


Selected Wide Flange Sections 





W SHAPES 
[6: p. 16] 
t, f 
Tr 
bi 
Elastic Compact 
Flange Web properties section 
Depth Width thickness Axis X-X criteria 

Designation d (in. by (in.) t, (In.) S (in.?) b;/2t, 
W 10X89 10.88 10.275 0.615 99.7 55 
W 10X60 10:25 10.075 0.415 67.1 7.38 
W10X49 10.00 10.000 0.340 54.6 8.96 
W10X45 10.12 8.022 0.350 49.1 6.49 
W 10X39 9.94 7.990 0.318 42.2 7.75 
W10X33 9.75 7.964 0.292 35.0 9.20 
W 10X25 10.08 5.762 0.252 26.5 6.70 
WI10X21 9.90 5.750 0.240 215 8.46 
W8X67 9.00 8.287 0.575 60.4 4.44 
W8X40 8.25 8.077 0.365 35.5 7.24 
WEX31 8.00 8.000 0.288 27.4 9.24 
W8X28 8.06 6.540 0.285 24.3 7.06 
W8X24 7.93 6.500 0.245 20.8 8.17 
W8X20 8.14 5.268 0.248 17.0 6.97 
W8X17 8.00 5.250 0.230 14.1 8.52 
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APPENDIX B 
EXAMPLE OF BEAM DESIGN PROCEDURE 


Parker [6: pp. 127-132] describes a general design procedure for a simple span, 
laterally supported beam with a uniformly distributed load. The beam is designed with 
A36 steel. The allowable deflection is limited to 1/360 of the span. Loading of the beam 


is uniformly distributed. 


Step 1: Compute the loads the beam will be required to support. For example, in a 
portion of a floor framing plan in Figure 1, a simple beam spans a distance 
of L = 10 feet and supports the floor a distance of 5 feet to either side of the 
beam. The floor load is 440 lb per sq ft, including the weight of the beam and 


other construction. 


Beam WA Column 





Figure 1 PORTION OF A FRAMING PLAN 


The uniformly distributed load per linear foot, w, including the weight on 


the beam, is: 
w = 5 ft x 440 lb per sq ft = 2,200 lb/ft = 2.2 kip/ft. 


The total uniformly distributed load, designated by W is the uniformly 
distributed load, w, times the span of the beam, L: 


WzwxL, 
or W = (2.2 kips/ ft) x (10 ft) = 22 kips. 


Step 2: Compute the reactions. Since the beam is syn  :trically loaded as shown in 
Figure 2., R, = R, = 22 kips/2 = 11 kips. 


W = 22 kips 





— 


‘111 kips 11 kips | Re 


rigure 2 REACTION DIAGRAM 


Ste J: Compute the maximum bending moment: 
M = WL/8 = (22 kips x 10 ft)/8 = 27.5 kip-ft. 
Step 4: Determine the required section modulus by use of the beam formula 
S = M/F,. where M is the bending moment in inch-pounds or kip-inches and 
Fyis the allowable bending stress in psi or ksi respectively. Since full 
lateral support is provided, the allowable bending stress of A36 steel is F, = 


24 ksi if a compact section is used. (ksi = kips per square inch). Then 


S = M/F, = (27.5 kip-ft x 12 in per ft)/24 ksi = 13.75 in. 
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Step D: 


Step 6: 


Step 7: 


Refer to the tables of structural shapes in Appendix A that give properties for 
designing, and select a beam with a section modulus equal to or greater than 
that required by Step 4. The beam selected will be adequate for bending 
stresses subject to checking of section classification (compact or 
noncompact) and L, or L, controls. Referring to Appendix A, [6: p. 16]. we see 
that a W8X17 has a section modulus of 14.1 in?, and we accept this as a trial 
beam. The compact section for the beam selection is checked. The compact 
section for A36 W-shapes is defined by AISC to be one in which the width- 
thickness ratio of the projecting compression flange, half-flange. does not 
exceed 52.2/VF,. The limiting value of the width-thickness ratio of the 
compression half-flange is computed to be 52.2/V36 = 8.70 for A36 steel 
[6: p. 22]. The compact section criteria in Appendix B shows that 
b,/2t, = 8.52 for the W8X17 section. Since this value is less than 8.7, the 
limiting value for A36 steel, the section is compact and the value of F, = 24 
ksi used in Step 4 is confirmed. 


Check the beam for shear by comparing the computed shearing stress 
f,= V/A, with the allowable F,. Checking the shear stress requires reference 
to Appendix A, which shows that d = 8 in and t,, = 0.23 in for the W8X17. 
Then since V = 11 kips: 


f, = V/A,= 11 kips/(8 in x 0.23 in) = 5.97 ksi, 
where A, = area of girder web in sq in. 


This value is less than the allowable 14.5 ksi, so the W8X17 is acceptable for 


shear. 


Compute the maximum deflection caused by the loading and compare it with 
the allowable delta, A. If the computed delta exceeds the allowable, the 
deflection formula may be solved for I required to limit delta to the 
allowable value [6: p. 83]. A beam section with an adequate moment of 
inertia may then be selected from the tables giving properties for designing. 
The allowable deflection is stated in the data as 1/360 of the span and is 
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(10 ft x 12 ft per in)/360 or 0.33 in. Since the allowable bending stress for 
this beam is 24 ksi, the following expression may be used to determine 


actual deflection [6: pp. 112-113, 131]: 


A= (5x f,x 2) / (48x Exo = (0.02483 L2) /d, 
A = (0.02483 x 102) / (8) = 0.31 in. 


The value is less than the allowable 0.33 in, so the W8X17 is acceptable for 
this loading. This expression is true only for simple steel beams with 
uniformly distributed loads. The extreme fiber stress, f,, is 24,000 psi for 
A36 steel, and E = 29,000,000 psi. 
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APPENDIX C 
NOMENCLATURE 
A partial selection of nomenclature and terms [6: pp. 2-5] 


American Institute of Steel Construction, Inc. 

Area of girder web (sq in) 

Modulus of elasticity of steel (29,000 kips per sq in) 

Bending stress permitted in the absence of axial force (ksi or psi) 

Allowable shear stress (ksi or psi) 

Specified minimum yield stress of the type of steel being used (ksi or psi). As 
used in AISC (American Institute of Steel Construction) Specification, “yield 
stress” denotes either the specified minimum yield point (for those steels that 
have yield point) or specified minimum yield strength (for those steels that do 
not have a yield point). 

Moment of inertia of a section (in?) 

Span length (ft) 

Maximum unbraced length of the compression flange at which the allowable 
bending stress may be taken at 0.66F, (ft) 

Maximum unbraced length of the compression flange at which the allowable 
bending stress may be taken at 0.60 F, (ft) 

Moment (kip-ft or kip-in) 

Reaction (kips) 

Maximum end reaction for 3-1/2 in. of bearing (kips) 

Elastic Section Modulus (in?) 

Statical (vertical) shear on beams (kips) 

Total] uniform load, including weight of beam (kips) 

Flange width of rolled beam or plate girder (in) 

Distance from the neutral axis to extreme fiber of beams (in) 

Depth of beam of plate girder (tn) 
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Computed bending stress (ksi or psi) 

Computed shear stress (ksi or psi) 

Actual unbraced length (in) 

Web thickness (in) 

Beam deflection (in) 

Uniformly distributed load per lineal foot (kip/lin ft) 


APPENDIX D 


ABBREVIATIONS 
(6: p. 5] 

Abbreviations Quantity 
cu ft cubic foot 
cu in cubic inch 
ft foot 
ft-lb foot-pound 
in inch 
in-lb inch-pounds 
kip 1000 pounds 
kip-ft kip-feet 
kip-in kip-inches 
ksf kips per square foot 
ksi kips per square inch 
lin ft linear foot 
Ib pounds 
lb per cu ft pounds per cubic foot 
lb per lin ft pounds per linear foot 
psf pounds per square foot 
psi pounds per square inch 
sq ft square foot 
sq in square inch 
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APPENDIX E 


MEASURES 


{17: pp. 120-126] 


/* measures.p created by hemant bhargava 1-20-89 */ 
/* Reference: Beyer's Mathematical Handbook */ 
/* measurement system, dimension, base unit, symbol, others */ 


/* Base units are fundamental or primary its. */ 
/* Units defined in terms of base units are tved units. */ 


/* base_unit(measurement_system, dimension, base_unit_symbol). */ 
base_unit(si,length,mt). 
base_unit(si,mass,kg). 
base_unit(si,time,sec). 
base_unit(si,current,amp). 
base_unit(si,temperature,degK). 
base unit(si, luminous intensity'.cd). 
base unit(si,' amount of substance',mol). 
base unit(si,currency,dollars). 
base unit(si,volume,.gal). 


base_unit(english, length. ft). 

base unit(english,mass.lb). 

base unit(english.,time,sec). 

base unit(english,current,amp). 
base_unit(english,temperature,degK). 
base_unit(english,'luminous intensity',cd). 
base unit(english,'amount of substance',mol). 
base unit(english,currency,pound). 


/* unit( measurement system, dimension, unit symbol). */ 
unit(si,length,mt). 
unit(si,length,cm). 
unit(si,length,km). 
unit(si,length,mm). 


unit(english ‚length ft). 

unit(english length in). 
unit(english,length,yd). 
unit(english,length,mi). 


unit(si,area,mt?2). 
unit(si,volume,mt43). 
unit(si,volume lt). 
unit(si,volume,gal). 
unit(si,volume,qt). 


unit(si,mass.kg). 
unit(si,mass,gm). 


unit(si,mass,mg). 
unit(si,mass,tons). 
unit(si,mass,ltons). 


unit(english,mass, lb). 
unit(english,mass.oz). 
unit(english,mass,stons). 


unit(_,time,sec). 
unit( .time,hr). 
unit( .time,min). 
unit( .time,day). 


unit(_.temperature.degK). 
unit(_.temperature,degC). 
unit(_,temperature,degF). 


unit(si,area,mt42). 
unit(si,volume,mt43). 
unit(si,volume.lt). 
unit(english,area.ft^2). 
unit(english ,volume.ft^3). 
unit(english,volume .floz). 


unit, code(si,2,mt). 
unit code(si,3,cm). 
unit code(si, 5, km). 
unit code(si, 7, mm). 


unit code(english. 1 1 ft). 

unit code(english, 13,1n). 
unit code(english.17 .yd). 
unit code(english, 19.m1). 


unit code(si.23. kg). 
unit code(si,29,gm). 
unit_code(si,31,mg). 
unit_code(si,37,tons). 
unit_code(si,41,ltons). 


unit_code(english,43,1b). 
unit_code(english,47,0z2). 
unit_code(english,53,stons). 


unit_code(_,61,sec). 
unit_code(_,67,hr). 
unit_code(_,71,min). 
unit_code(_.73,day). 


unit_code(_.83.degK). 
unit_code(_,89,degC). 
unit code( .91.degF). 
unit_code(_,101.dollars). 


/* The convention is to state across-system laws 
only for the base unit in each system 
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e.g. between mt and fî ut not cm and ft. 
within system conversion is done using the conversion 
laws with the base unit as the meeting point */ 


/* Base conversion laws (ACROSS systems) */ 
base_conv_law(mt = 0.3048 * ft). 
base_conv_law(kg = 0.45359237 * lb). 


/* RHR */ 
/* Conversion from base unit of a system to other unit in it */ 
/* The query we are interested in answering is 
C=?F =? in: conv_law(Ul =C +F * U2). */ 
/* length */ 
conv_law(mt = 100 * cm). 
conv_law{mt = 0.001 * km). 
/* conv_law(mt = 1609.344 * mi). */ 
conv law(mt = 1000 * mm). 


/* length */ 

conv_law(ft = (1/3) * yd). 
conv_law(ft = (1/5280) * mi). 
conv law(ft 2 12 * in). 


/* mass */ 

conv, law(kg = 1000 * gm). 

conv law(kg 2 1000000 * mg). 
conv law(kg = (1/1000) * tons). 
conv law(kg 2 (1/1000) * Itons). 


/* mass */ 
conv. law(lb = 16 * oz). 
conv. law(Ib 2 (1/2000) * stons). 


/* time */ 

conv_law(sec = (1/60) * min). 
conv law(sec = (1/3600) * hr). 
conv_law(sec = (1/86400) * day). 


/* temperature */ 

conv_law{degK = 273.15 + degC). 
conv_taw(degF = 32 + (9/5) * degC). 

conv. law(degC = (-32 * 5 / 9) + (5/9) * degF). 


conv_law(mt 3 = 1000 * It). 
conv_law(gal = 4 * qt). 


conv_law(U = 0 + 1 * U). 
conv_law(U14A = O + Factor * (U2*A)) :- 
conv_law(Ul = Factor * U2), 
Factor is pow(F,A). 


/* 969696969696969690690696909096909/0906906969090969696909069696969696969696969069096909696 */ 
/* For converting between base unit */ 


conv_law(BU1 = 0 + F12 * BU2) :- 
base_conv_law(BU1 = F12 * BU2). 
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conv_law(BU1 = 0 + F12 * BU2) :- 
base conv law(BU2 = F21 * BUI). 
]12 15 1 / "21 


conv_law(BUl = Fl2 + 1 * BU2):- 
base_conv_law(BU1 = F12 + BU2). 
conv. law(BU1 z F12 € 1 * BU2) :- 
base conv law(BU2 2 F21 « BUI), 
F12 is (-1) * F21. 


/* Within system: base unit to other unit*/ 
conv_law(BU = 0 + F * U) :- 

conv_law(BU = F * U). 
conv_law(BU = C + 1 * U):- 

conv law(BU = C + U). 


/* For arbitrary conversions (within or across systems) */ 
conv_law(Ul = 0 + F12 * U2) :- 
unit(S1,Dimension,U 1), 
base unit(S1.Dimension,BU 1), 
conv_law(BU1 =0+F1 * Ul), 
unit(S2,Dimension, U2), 
base_unit(S2,Dimension,BU2), 
conv_law(BU2 = O + F2 * U2), 
conv law(BUl = 0 + BF12 * BU2), 
F12 is F2 * BF12 / F1. 


/* For arbitrary conversions (within or across systems) */ 
conv_law(U1 = F12 + 1 * U2) :- 
unit(S 1, Dimension,U 1), 
base unit(S1,Dimension,BU 1), 
conv_law(BUl] =F1+1* Ul), 
unit(S2,Dimension, U2), 
base_unit(S2,Dimension,_,BU2), 
conv_law(BU2 = F2 + 1 * U2), 
conv_law{BU1 = BF12 + 1 * BU2), 
F12 is F2 + BF12 - Fl. 


unit_alias(meters,mt). 

unit alias(kilograms.kg). 
unit_alias(seconds, sec). 
unit_alias(amperes, amp). 
unit_alias('degree Kelvin',degK). 
unit_alias(candela,cd). 
unit_alias(mole,mol). 
unit_alias(dollars,dollars). 


/* convert from units 2 to units 1 */ 
convert(Number,Units1,Number,Units1) :- !. 
convert(Number2,Units2, Number] ,Units1) :- 

conv law(Units2 - F1 + F2 * Units1), 
Number] is F1 « (F2 * Number2). 
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LAWS FOR DIMENSIONAL CONSISTENCT 


Fe ame 
249 
- - -e و - - جي‎ 
2 È ie m ma - — æ - - سے صم مر‎ = 
== “<= Ze Aa a = CERES HK Tw SR ru a META C Wd. 2 Caf SE RE" S = © ESE IE: c 
u - س سس و‎ > 5 — © © DA wann ® ao چ ج که‎ ee Poe ee — — T © è è tè —. = pe —--p  è& al nn 
= * © < ^ - + 
z omm re as UE = -— LES e ==” 3 irat uw X 
— m ee - -— = ans a a. ——— —— A - = - م‎ 
x 2 - 
rm sf lé Lees E 2 
ess +" se dda سد‎ o 
1-527 & BE "1 v2 ممه‎ o z x 
e ee w = 9 ن د‎ A— A-0009 j 
کڪ‎ T 2089 4 Ww ER De c acm 9 = NF Y è i csi cali Di: E uw ^ -—— — = PE o 3 sos - 9 = 
t + سر‎ * e — e "-—-] ^? StL. = > 
nt [c— —— © = - a - mit ur ts 2e ^ o uw 9 o ais ow AA مهمه‎ - oo =. —— - Nur عياب نت‎ y os o 
= 
con ulm YN a sr 4 nmel & yd 
ھم سے‎ A. 9 - «= — e RITI BD - æ دم‎ LE 
sms frs pe: ye 
< a سے‎ tan ad o u o 
~ a, 
17-77 Dmi ) a ml = 2m 
ے ست‎ "7 9 ^ = a uf Ce ur 
- 
dim validi 2-5) -- rite : 
هم حسم ہے‎ -— T e `. æ a © = i. — en. ^ - - > 
> if af e mw ox* 
cM NU NE. mt © LA o 
— 
LEFT = Es ) 4 es m= - = 
= ee ee è on — «è رس‎ > 
“=. - t -— À - Le a 
dim velid[M. AB) :- dim velid(M. : 
—— 9 AL. © © 7 e ———— oa P... ES 
2 € و« حم‎ 
dim valid(M =) 
1 ica ao u | 
ze Le z * = sa» 
i men qo i 
MY di ra rl Né ~ 
مه مص‎ e ¿rr A 0 de 


à. 1: ina à A È Sie = * ¢ % 
im_valid(M. A*B) :- dim_ alid(M. A). 
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dim vaitd(M. 


dim valád(M.sumi . X)) :- dim valid(M. X). 


dim valid(M.prod( . X)) :- dim_valid(M. XI. 
dim valid(M. AzB) :- units{M. A, DA), 
units(M, B, DB). 
units eq(DA. DB, true). 


cim valid(M, A<B) :- units(M, A. DA). 
units{M, B, DB). 
units_eq(DA. DB. true). 


dim_valid(M. A>B) :- units(M. A. DA). 
units(M. B. DB). 
units eq(DA. DB, true). 


dim, valid(M. A=<B) > units!" 
units(M. B. DB). 


A. DA). 


Y 


units_eq(DA, DB, true). 
dim_valid(M, A>=B) :- units(M, A, DA), 
units(M, B, DB), 
units_eq(DA, DB, true). 
dim_valid(M, Exp, false). 
/* Laws for Computing derived units */ 
units(M, A+B, UA) :- units(M, A, UA), 
units(M, B, UB), 
units_eq(UA, UB, true). 
units(M, A-B, UA) :- units(M, A, UA), 
units(M, B, UB), 
units_eq(UA, UB, true). 
units(M, A*B, U) :- units(M, A, UA), 
units(M, B, UB), 
simplify(UA*UB, U). 
units(M, A/B, U) :- units(M, A, UA), 
units(M, B, UB). 
simplify(UA/UB, U). 
units(M, A^B, U) :- units(M, A, UA), 
units(M, B, UB), 
simplify(UA^UB, U). 
units(M, pow(A, B), U) :- units(M, A, UA), 
units(M, B, UB). 
simplify(UA^UB, U). 
units(M ,max(_,X), U) :- units(M, X, U). 
units(M,min(_,X), U) :- units(M, X, U). 
units(M,sum( ,X), U) :- units(M, X, U). 
/* units(M, prod(_,X), U) :- units(M, X, U). requires multiplication of each */ 
units(M, Local. var, U) :- local varí(. , M, Local var, ,U). 
units eq(A, A, true). 
/* If any of the definitional expressions is dimensionally invalid, the model is not d.c. */ 


/* If a model is d.c. the second argument of the dim valid model predicate will be 'true'. 
If it is not, the second argument will contain the invalid expression. */ 


dim_valid_model(Model, DefE) :- defE(Model, _, DefE), 
not(dim_valid(M, DefE)), !. 


/* If any of the conditional expressions is dimensionally invalid, the model is not d.c. */ 
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dim, valid. model(Model, ConE) :- defE(Model, , ConE). 
not(dim. valid(M, ConE)), !. 


/* Else it must be OK */ 


dim, valid. model(Model, true). 


APPENDIX G 
A NUMERIC METHOD 


/* numeric.p */ 

/* 12 Feb 1991 */ 

/* Get numeric value for dertved units. */ 

/* Assume database of prime values for fundamental units*/ 


/* numeric(Dimensional expression, Numeric value) */ 


numeric(Exp l*'Exp 2, Numeric value) :- 
numeric(Exp 1. N1), 
numeric(Exp 2, N2), 
Numeric value is N1*N2. 


numeric(Exp 1/Exp 2, Numeric value) :- 
numeric(Exp 1, N1), 
numeric(Exp 2, N2), 
Numeric value is N1/N2. ١ 


numeric(Exp_1+Exp_2, Numeric_value) :- 
numeric(Exp 1. N1), 
numeric(Exp 2, N2), 
N1 = N2, 
Numeric_value is N1. 


numeric(Exp 1-Exp 2, Numeric value) :- 
numeric(Exp 1, N1), 
numeric(Exp 2. N2), 
N1 = N32, 
Numeric value is Ni. 


numeric( Exp 1^Exp 2, Numeric value) :- 
numeric(Exp 1, N1), 
numeric(Exp. 2, N2), 
N2 =0, 
Numeric_value is Nl. 
numeric(Unit, Prime) :- unit_code(System, Prime, Unit). 


numeric(X, 1) :- number(X). 
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APPENDIX H 


A LIST PROCESS 


/* multiply list([[A]! [B]], Answer). */ 
multiply list([]. 1). 


multiply list([Al B], Answer) :- 
multiply list(A, X). 
multiply list(B, Y). 
Ar "eris XY. 


multiply A, Ans ) :- 
nu:..»er(A), 
Ans is A. 


/* divide list([[A] | [B]], Answer). */ 
divide list((]. 1). 


divide list([Al B], Answer) :- 
divide list(A, X), 
divide list(B, Y), 
Answer is X/Y. 


divide list(A, Ans ) :- 
number(A), 
Ans is A. 


/* add list */ 

/* add list([[A] ! [B. C, D]]. Answer). */ 
/* add list([Al [B, C, D]]. Answer). */ 
add list([], O). 


add list([Al B], Answer) :- 
add list(A, X). 
add list(B, Y), 
Answer is X+Y. 


add_list(A, Answer) :- 
number(A), 
Answer is A. 


/* subtract list([[A] | [B]], Answer). */ 
subtract list([], O). 


subtract list([Al B], Answer) :- 
subtract list(A, X), 
subtract list(B, Y). 
Answer is X-Y. 


subtract_list(A, Ans ) :- 
number(A), 
Ans is A. 
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APPENDIX I 


BEAM MODEL 


/* model(model_name, description). */ 
/* model_source(model_name, r >l_reference). */ 
/* variable(variable_name, desc on). */ 
/* var quiddity(variable name, a of quantity, quiddity, storage units). */ 
/* defE(model_nar =, expression_number, definitional expression). */ 
/* local. var(model name, variable r :23e, local. name, «exo or endo», units). */ 
/* setDef[model name, index name,  escription', index range). */ 
/* indexvar(index, model name. 'descnption', [type(index), index rangel). */ 
/* in. model(model name, use(model name, local variable), 
to eval(local variable)) */ 
/* run report(model name, [variable list]). */ 
/* scenario(model name(scenario number), description). */ 
/* m, s pair(model name, model name(scenario number)). */ 
/* datum(variable name, model name(scenario number), data). */ 


/* Variable base */ 


var quiddity(allowable shear. stress, mass/area, *, kip/1n^2). 
variable(allowable, shear. stress, 'Allowable shear stress of beam’). 


var quiddity(actual deflection, length, *, in). 
variable(actual deflection, 'Actual deflection of beam’). 


var_quiddity(allowable_«.. ‘lection, length, *, in). 
variable(allowable_deflection, 'Allc:’ able deflection of beam’). 


var_quiddity(beam_depth#i, length, *, in 
var(beam_depth#i, ['Depth of beam , il). 


var quiddity(bending stress, mass/area, *, kip/in^2). 
variable(bending stress, 'Bending stress). 


var quiddity(computed, section modulus, volume, *, in^3). 
variable(computed section modulus, 'Computed section modulus from beam 
formula S=M/Fb)). 


var_quiddity(computed_shear_stress, mass/area, *, kip/in?2). 
variable(computed_shear_stress, 'Computed shear stress’). 


var_quiddity(depth_solution, length, *, in). 
variable(depth_solution, 'Depth of wide flange solution’). 


var_quiddity(length_of_beam, length, *, ft). 
variable(length_of_beam, 'Length of beam’). 
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var_quiddity(max_vertical_shear, mass, *, kip). 
variable(max_vertical_shear, 'Maximum vertical shear of beam’). 


var quiddity(moment of beam, mass*length, *. kip"ft). 
variable(moment of beam, 'Moment of beam’). 


var_quiddity(reaction, mass, *, kip). 
variable(reaction, ‘Reaction of beam’). 


var quiddity(section modulus?i, volume, *, in^3). 
variable(section_modulus#i, [Table section modulus ', 1)). 


var_quiddity(uniformly_distributed_load, mass/length, *, kip/ft). 
variable(uniformly_distributed_load, ‘Uniformly distributed load’). 


var_quiddity(w_shape#i, _, *, _). 
variable(w_shape#i, [Wide flange shape ', il). 


var_quiddity(web_thickness#i, length,*, in). 
variable(web_thickness#i, ['Web thickness of beam ', iJ). 


var quiddity(web thickness solution, length,*, in). 
variable(web thickness solution, "Web thickness of wide flange solution’). 


var_quiddity(weight_of_load, mass, *, kip). 
variable(weight_of_load, ‘Weight of load’). 


var_quiddity(wide_flange_solution, _, *, _). 
variable(wide flange solution, 'Wide flange solution"). 


model(beam, 'beam: Model for selection of simple span beams. ). 
model_source(beam,'Parker: Simplified Design of Structural Steeel.'). 
m_s_pair(beam, beam(1)). 


/* step_1: Compute weight */ 
model(step_1, 'Compute weight of floor load.’). 
model_source(step_1,'Parker'). 
m s pair(step 1, beamí(1)). 

local var(weight of load, step 1, cap_W. endo, kip). 

local var(length of beam, step 1. cap L. exo, ft). 

local var(uniformly distributed load, step 1, w, exo. kip/ft). 
/* weight of load(kip) 2» uniformly distributed load(kip/ft) * length of beamíft) */ 
defE(step_1, 1, cap_W =: w * cap_L). 

in_model(step_1, use(step_1, cap_W), to_eval(cap_W)). 
run, report(step 1, [cap W, w, cap Ll). 


41 


/* step_2: Compute reactions */ 

model(step_2, ‘Compute beam reactions. ). 

model_source(step_2,'Parker'). 

m_s_pairístep_2, beam(1)). 
local_var(reaction, step_2, cap_R, endo, kip). 
local_var(max_vertical_shear, step_2, cap_V, endo, kip). 
local_var(weight_of_load, step_2, cap_W, exo, kip). 

/* reaction(kip) = (weight_of_load(kip) / 2) */ 

defE(step_2, 2-1, cap_R =: (cap_W / 2)). 

/* max_vertical_shear(kip) = reaction(kip) */ 

defE(step 2. 2-2, cap V =: cap_R). 
in model(step 2. use(step 1, cap W), to eval(cap W)). 
in model(step 2, use(step 2, cap RJ), to eval(cap. R)). 
in model(step 2. use(step 2, cap V), to eval(cap. V)). 

run report(step 2, [cap R., cap. Vl). 


/* step 3: Compute moment */ 
model(step 3, ‘Compute moment.’). 
model source(step 3,' Parker). 
m s pair(step 3, beam(1)). 
local var(moment of beam, step 3, cap M, endo, kip*ft). 
local var(length of beam, step 3, cap L, exo, ft). 
local var(weight of load, step 3, cap W, exo, kip). 
/* moment of beam(kip*ft) 2 (weight of load(kip)) * length of beam(ft))/ 8 */ 
defE(step 3, 3, cap M z: (cap W * cap L)/8). 
in model(step 3, use(step 1, cap W), to eval(cap. W)). 
in model(step. 3, use(step 2, cap R). to eval(cap. R)). 
in model(step 3, use(step 2, cap V), to eval(cap. V)). 
in model(step 3, use(step 3, cap M). to eval(cap. M)). 
run report(step 3, (cap MJ). 


/* step 4: Compute section modulus */ 
model(step 4, 'Compute section modulus). 
model source(step 4,' Parker). 
m s palir(step 4, beam(1)). 
local var(computed, section modulus, step 4. cap. S, endo, in^3). 
local var(moment of beam, step 4. cap M, exo, kip*ft). 
local var(bending stress, step 4, cap Fb, exo, kip/in^2). 
/* computed section, modulus(in^3) 2 (moment of beam(kip*ft) * (12 
in/ft))/bending stress(kip/in^2)) */ 
defE(step 4, 4, cap S z: (cap M *12) / cap Fb). 
in model(step 4, use(step 1, cap W), to eval(cap W)). 
in model(step 4, use(step 2. cap RJ), to eval(cap. R)). 
in model(step 4, use(step 2, cap V). to eval(cap. V)). 
in model(step 4, use(step 3. cap M). to eval(cap. M)). 
in model(step 4, use(step 4, cap. S). to eval(cap. S)). 
run report(step 4, [cap S. cap Fb]). 
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/* step_5: Select beam */ 

model(step_5, ‘Select beam using computed section modulus). 

model. source(step. 5.'Parker). 

m, s pair(step. 5, beam(1)). 

06 local var(solution list, step 5, [list], endo, ). 
local var(test variable, step. 5, test, endo, ). 
local, var(wide flange solution, step 5. w sol, endo, _). 
local var(depth solution, step. 5, d, endo, in). 
local var(web thickness, solution, step 5, tw, endo, in). 
local, var(computed. section modulus, step 5. cap. S, exo, in^3). 
local var(section modulus, step. 5, s mod. exo, in^3). 
local var(w. shape?sti, step 5, shapetti, exo, _). 
local var(beam depth?fti, step 5, depth#i, exo, in). 
local var(web thickness?i, step 5, thickfsi, exo, in). 

9o defE(step. 5, 5-1, listz: vector(i in a, list, 1 «-- s modfti »- cap. S)). 

96 defE(step. 5, 5-2. w sol =: min_list(list)). 

defE(step. 5, 5-3, w sol z: shapes 15). 

defE(step_5, 5-4, d =: depth#15). 

defE(step_5, 5-5, tw=: thick#15). 

% setDef(step_5, a_list, 'Possible beam choices', 1...15). 

% indexvar{i, step_5, ‘Beam choice index‘, [i in a_list]). 
indexvar{i, step_5, ‘Beam choice index’, {integer{i), 1... 15)). 
in_model(step_5, use(step_4, cap_S), to_eval(cap_S)). 
in_model(step_5, use(step_5, w_sol), to_eval(w_sol)). 
in_model(step_5, use(step_5, d), to_eval(d)). 
in_model(step_5, use(step_5, tw), to_eval(tw)). 

% in_model(step_5, use(step_5, list), to_eval(list)). 
in_model(step_5, use(step_5, test), to_eval(test)). 

run_reportístep_5, [w_sol, d, tw]). 


/* Step_6: Check shear stress */ 
model(step_6, 'Check shear stress”. 
model_source(step_6,'Parker'). 
m s pair(step 6, beam(1)). 
local var(computed shear. stress, step 6, fv, endo, kip/in^2). 
local var(shear check, step 6, shear, endo, _). 
local var(allowable shear stress, step 6. cap Fv. exo, kip/in^2). 
local var(max vertical shear, step. 6, cap. V, exo, kip). 
local var(depth solution, step 6. d, exo, in). 
local var(web thickness solution, step 6, tw, exo, in). 
/*computed, shear. stress(kip/in^2)2max vertical shear(kip)/(depth(in) * thick(in)).*/ 
defE(step 6, 6-1, fv z: cap V / (d * tw)). 
/* If allowable shear. stress(kip/in^2) »2- computed shear. stress(kip/in^2) then 
(shear. check = ok) */ 
defE(step_6, 6-2, shear =: (ok <-- cap_Fv >= fv)). 
in_model(step_6, use(step_2, cap_V), to_eval(cap_V)). 
in_model(step_6, use(step_5, d), to_eval(d)). 
in_model(step_6, use(step_5, tw), to_eval(tw)). 
in_model(step_6, use(step_6, fv), to_eval(fv)). 
in_model(step_6, use(step_6, shear), to_eval(shear)). 
run_report(step_6, [shear, fv, cap Fvl). 


/* Step_7: Check deflection */ 
model(step_7, 'Check deflection’). 
model_source(step_7,'Parker'). 
m s pair(step 7, beam(1)). 
local var(allowable deflection, step 7, allow, endo, in). 
local var(actual deflection, step 7, actual, endo, in). 
local var(deflection check, step 7, defl, endo, ). 
local var(length of beam, step 7, cap L, exo. ft). 
local var(depth solution, step. 7, d, exo, in). 
/* allowable deflection(in) 2 (length of beamt(ft) * (12 in/ft))/ 360) */ 
defE(step. 7, 7-1, allow z: (cap. L * 12) / 360). 
/* actual deflection(in) 2 (0.02483(1in^2/ft^2) * length of beam(ft^2) / d(tn) */ 
defE(step. 7, 7-2, actual =: (0.02483 * cap L^2) / d). 
/* If actual deflection(in) 2« allowable. deflection(in) then (deflection check 2 ok) */ 
defE(step_7, 7-3, defl =: (ok <-- actual =< allow)). 
in_model(step_7, use(step_5, d), to_eval(d)). 
in_model(step_7, use(step_7, allow), to_eval(allow)). 
in_model(step_7, use(step_7, actual), to_eval(actual)). 
in_model(step_7, use(step_7, defl), to_eval(defl)). 
run report(step 7, [defl, actual, allow]). 


scenario(beam(1), 'This is a test scenario of the beam model). 


datum(allowable, shear. stress, beam(1), 14.5). 
datum(bending, stress, beam(1), 24.0). 
datum(length of beam, beam(1). 10.0). 
datum(uniformly distributed load, beam(1), 2.2). 


/* Table 1-1, Parker, Simplified Design of Structural Steel */ 


datum(w. shapes 1, beam(1), w10 89). 
datum(w shapef2, beam(1), w10 60). 
datum(w. shapes3, beam(1), w10. 49). 
datum(w_shape#4, beam(1), w10_45). 
datum(w_shape#5, beam(1), w10_39). 
datum(w_shape#6, beam(1), w10_33). 
datum(w_shape#7, beam(1), w10_25). 
datum(w_shape#8, beam(1), w10_21). 
datum(w_shape#9, beam(1), w8_67). 

datum(w_shape#10, beam(1), w8_40). 
datum(w_shape#11, beam(1), w8_31). 
datum(w_shape#12, beam(1), w8_28). 
datum(w_shape#13, beam(1), w8_24). 
datum(w_shape#14, beam(1), w8_20). 
datum(w_shape#15, beam(1), w8_17). 


datum(section_modulus#1, beam(1), 99.7). 
datum(section_modulus#2, beam(1), 67.1). 
datum(section_modulus#3, beam(1), 54.6). 
datum(section_modulus#4, beam(1), 49.1). 
datum(section_modulus#5, beam(1), 42.2). 
datum(section_modulus#6, beam(1), 35.0). 
datum (section_modulus#7, beam(1), 26.5). 


datum(section_modulus#8, beam(1), 21.5). 
datum(section_modulus#9, beam(1), 60.4). 


datum(section_modulus# 10, beam(1), 35.5). 
datum(section_modulus#11, beam(1), 27.4). 
datum(section_modulus#12, beam(1), 24.3). 
datum(section_modulus#13, beam(1), 20.8). 
datum(section_modulus#14, beam(1). 17.0). 
datum(section_modulus#15, beam(1), 14.1). 


datum(beam_depth#1, beam(1), 10.88). 
datum(beam_depth#2, beam(1), 10.25). 
datum(beam_depth#3, beam(1), 10.00). 
datum(beam_depth#4, beam(1), 10.12). 
datum(beam_depth#5, beam(1), 9.94). 

datum(beam depth£6, beam(1), 9.75). 

datum(beam depth£27, beam(1), 10.08). 
datum(beam, depth£8, beam(1), 9.90). 

datum(beam depth£49, beam(1), 9.00). 

datum(beam depth? 10, beam(1), 8.25). 
datum(beam depth?& 11, beam(1), 8.00). 
datum(beam, depth? 12, beam(1), 8.06). 
datum(beam, depth? 13, beam(1). 7.93). 
datum(beam, depth? 14, beam(1), 8.14). 
datum(beam depth£& 15, beam(1). 8.00). 


datum(web. thickness? 1, beam(1), 
datum(web_thickness#2, beam(1), 
datum(web_thickness#3, beam(1), 
datum(web_thickness#4, beam(1), 
datum(web_thickness#5, beam(1), 
datum(web_thickness#6, beam(1), 
datum(web_thickness#7, beam(1), 
datum(web_thickness#8, beam(1), 
datum(web_thickness#9, beam(1), 


0.615). 
0.415). 
0.340). 
0.350). 
0.318). 
0.292). 
0.252). 
0.240). 
0.575). 
datum(web_thickness#10, beam(1), 0.365). 
datum(web. thickness? 11, beam(1), 0.288). 
datum(web. thickness? 12, beam(1), 0.285). 
datum(web thickness? 13, beam(1), 0.245). 
datum(web thickness? 14, beam(1), 0.248). 
datum(web thickness? 15, beam(1), 0.230). 
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APPENDIX J 


SAMPLE OUTPUT 


COMPUTE WEIGHT 
TEFA»»run(step. 1, beam(1)). 


cap W222.0 kip w22.200000 kip/ft cap L-10.0 ft 


COMPUTE REACTIONS 
TE "A»»run(step 2. beamí1)). 


cap_R=11.0 kip cap_V=11.0 kip 


COMPUTE MOMENT 
TEFA>>run(step_3, beam(1)). 


cap_M=27.500000 kip*ft 


COMPUTE SECTION MODULUS 
TEFA»»run(step 4, beam(1)). 


cap S213.750000 in^3 cap Fbz24.0 kip/in^2 


SELECT BEAM 
TEFA»»run(step. 5, beam(1)). 


w_sol=w8_17 d=8.0in tw=0.230000 in 


CHECK SHEAR 
TEFA>>run(step_6, beam(1)). 


shear=ok fv=5.978260 kip/in‘2 cap_Fv=14.500000 kip/in?2 


CHECK DEFLECTION 
TEFA>>run(step_7, beam(1)). 


defl=ok actual=0.310375 in allow=0.333333 in 
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APPENDIX K 


DESCRIPTION OF COMPUTER PLATFORM 


A. SYSTEM DEVELOPMENT INFORMATION 
1. Computer platform: Macintosh IIsi. 
2. Operating system: 6.07. 
3. Processor: Motorola 68030, 20 mHz. 
4. RAM: 5 megabytes. 
5. Hard disk storage: 80 megabytes. 


B. SYSTEM REQUIREMENTS 
1. Computer platform: Mac II series. 
2. Operating system: 6.05 or later. 
3. RAM: 4 megabytes, minimum. 
4. Hard disk: Highly desirable. 


C. PROGRAMMING ENVIRONMENT 
1. Model Management System: TEFA, version 1.10. 
2. Programming language: Advanced A.I. Systems Prolog 2.0g. Prolog is the 
embedded language used by TEFA. 
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